AD-A147  435 


REGION  1  BIRKELRND  CURRENTS  FROM  A  GLOBAL  MHD 
SIMULATION  OF  THE  MAGNETOSPHERES)  NAVAL  RESEARCH  LAB 
WASHINGTON  DC  J  A  FEDDER  ET  AL.  28  SEP  84  NRL-MR-5297 

F/G  4/1 


1/1 


UNCLASSIFIED 


NL 


This  work  was  supported  py  tl 
and  the  National  Aeronautics 


ft 


ff 


m 


$ 


j.  A.  fedder  and  J.  g.  Lyon 

Geophysical  and  Plasma  Dynamics  Branch 
Plasma  Physics  Division 


NAVAL  RESEARCH  LABORATORY 


Wasfctastan,  D.C 


AmmwsA  tor  pubic  release;  Orttibrnta  nflaML 


11  07  04  fr* 


REPORT  DOCUMENTATION  PAGE 


UNCLASSIFIED 


:«V\G  DRGANZA'ON  RERCR*  N_V3ER;S; 
NRL  Memorandum  Report  5397 


x  aoc5HSS  Cry  Sure  and  ZIP  Code) 
Washington,  DC  20375 


5  MON  "CR-NG  DRGANiZA-fN  REPORT  NUM3ERIS) 


MCNiTDRi.NG  ORGANIZATION 


*b  ADDRESS  Cry,  Stare,  and  ZIP  Cade) 


3a  NAVE  QF  =GNO'NG  ,  SPONSORING 
ORGAVZA*  ON 


3o  OP=iC£  SVM80L 
i if  applicable) 


9  PROCUREMENT  NSTRUN1ENT  OENTlf iCATION  NUMBER 


ONR  and  NASA 


3<  ADDRESS  (C/fy,  State  and  ZIP  Code) 

Arlington,  VA  22217  Greenbelt,  MD  20771 


•o  SOURCE  OF  cUNDlNG  NUM6ERS 


PROGRAM 
ELEVEN"  NO 

61153N 


PROJECT 

task 

NO 

*0  RR033- 

W-15494 

02-44 

WORK  UNIT 


"■  -E  include  security  Clarification) 

Region  1  Birkeland  Currents  from  a  Global  MHD  Simulation  of  the  Magnetosphere 


PERSONAL  AvTHORiS) 

Fedder,  J.A.  and  Lyon,  J.G. 


'3a  "V?E  OF  REPORT 

Interim 


13b  "IME  COVERED  1 A  DATE  OF  REPORT  lYear,  Month.  Day) 

prom  10/84  ro  10/85  1984  September  28 


'6  Supplementary  no’jton  This  work  was  supported  by  the  Office  of  Naval  Research  and  the 
National  Aeronautics  and  Space  Administration. 


cosati  cooes 


GROUP 


’8  SUBJECT  "E.RMS  Continue  on  reverse  it  necessary  and  identify  by  block  number) 

MHD  simulation  Birkeland  current 

Magnetosphere  Polar  current  systems 


'3  A3S’RAC~  Continue  on  reverse  ,f  necessary  and  identify  by  block  number) 

Global  numerical  MHD  simulation  data  for  the  earth’s  magnetosphere  has  been  used  to  calculate 
Birkeland  current  systems  in  the  polar  ionosphere.  In  this  report,  we  show  the  results  for  the  Region  1 
Birkeland  currents.  We  demonstrate  how  the  currents  are  organized  on  open  or  closed  field  lines  in 
daytime  and  nighttime  sectors.  The  current  source  regions  are  identified  and  driving  mechanisms  are 
discussed. 


30  O'VR'SU'ON  AvAILASlL'Y  Or  A3STSAC" 

S  ^NC-ASS  F  ED  DNL  M1TED  □  SAVE  AS  RPT 


Z2»  NAME  OF  RESPONSBlS  NDiV'DUAL 

•J.  A.  Fedder 


2’  abstract  security  classification 
□  DTic  USERS  UNCLASSIFIED 


225  TE-E?uONt  (include  Area  Code)  22c  CPF'Ct  >YM0OL 

(202)  767-3630  Code  4780 


10.  SOURCE  OF  FUNDING  NUMBERS 


WORK  UNIT 
ACCESSION  NO. 

DN330-1S9 

DN3S0-475 


7^ 


,*J 


:-j 


CONTENTS 


INTRODUCTION  .  1 

METHOD  .  2 

RESULTS  .  4 

DISCUSSION  AND  CONCLUSIONS  .  5 

ACKNOWLEDGMENTS .  9 

REFERENCES  .  14 


c 

M 


a 


DTIC 

ELECTEf* 
NOV  1  4  1984  ;  •  1 

B 


in 


Dl3t 


M 


i  Access i 

V 

!  NT « S  r 
1  DTIC  I. 

Ur.  .'•a-. 

:  Ju-n1'* 


Hy. 

Dinti i 
A  w '  ’ 


/■ 


-1 


1 


4 


vlvIy-VliU 


•y-vCz-vS 


VV 


REGION  1  BIRKELAND  CURRENTS  FROM  A  GLOBAL 
MHD  SIMULATION  OF  THE  MAGNETOSPHERE 


Introduction 

-  - -T'The  study  of  large  scale  geomagnetic  field  aligned  currents,  which 

couple  the  dynamics  of  the  magnetosphere  and  the  high  latitude  ionosphere, 
has  been  an  area  of  intense  research  interest  in  recent  years.  These 
currents,  first  postulated  by  Birkeland  [1908]  to  be  responsible  for 

geomagnetic  activity  and  auroral  phenomena  in  the  polar  regions.,  were  first 
observed, by  Arnoldy^  [1974],  Zmuda  and  Armstrong  [1974a, b],/  and  Sugiura 
[1975T  and  have  now  been  extensively  measured  and  documented  by  satellite 
and  rocket  magnetometers .  (  Much  of  the  recent  experimental  research  has 

been  reviewed  by  Potemra  [1979],  Potemra  et  al.  [1979],  Burch  and  Heelis 

[1978],  Saflekos  et  al.  [1982],  Greenwald  [1982],  Burke  [1982],  and  Potemra 
[1982]. 

The  large  scale  Birkeland  currents  are  observed  in  two  rings  which 

encircle  the  polar  regions.  The  poleward  ring  has  been  called  the  "Region 
1"  system  and  the  equatorward  ring  the  "Region  2"  system,,  [Iijima  and 

Potemra,  1976a, b].  The  Region  I  currents  are  earthward  in  the  morning 
sector  and  away  from  the  earth  in  the  evening  sector,  while  the  Region  2 
currents  are  in  the  opposite  direction  in  each  sector.  The  Region  1 
currents  have  been  identified  as  magnetospher ic  driving  currents  and  are 

directly  associated  with  the  solar  wind-magnetospher ic  dynamo,  whereas  the 
Region  2  currents  are  considered  to  be  response  currents  [Schield  et  al. , 

1969;  Wolf,  1974].  The  Region  1  and  Region  2  Birkeland  current  systems  and 
their  sources  and  sinks  are  a  fundamental  link  in  the  solar  wind- 
magnetosphere-ionosphere  interaction  and  dynamics. 

Manuscript  approved  June  6,  1984. 
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Quantitative  theoretical  and  simulation  studies  of  Birkeland  current 
systems  and  magnetospheric  dynamics  were  begun  by  Wolf  and  Harel  [1979]  and 
Harel  et  al.  [1979].  This  pioneering  work  and  continuing  research  [Harel 
et  al.,  1981a;  Harel  et  al.,  1981b;  Spiro  et  al.,  1981;  Wolf  et  al.,  1982; 
Chen  et  al.,  1982;  and  Karty  et  al.,  1982]  has  led  to  a  major  advance  in 
understanding  Birkeland  currents  and  their  role  in  magnetospheric 

dynamics.  These  simulation  studies  have  treated  the  inner  magnetosphere- 

ionosphere  system  in  detail  and  therefore  have  provided  a  quantitative 
model  for  the  Region  2  currents  which  are  consistent  with  observations. 
However,  the  Region  1  currents  are  outside  or  on  the  boundary  of  the 
simulated  region  of  the  magnetosphere  and  are  therefore  untreated. 

In  this  brief  report  we  present  preliminary  simulation  results  for  the 
Region  1  Birkeland  current  systems.  The  results  are  extracted  from  a 
global  MHD  solar  wind  magnetosphere  simulation  and  therefore  provide  a 
description  of  the  driving  current  system.  Owing  to  limited  numerical 
resolution  and  an  inadequate  treatment  of  the  ionosphere,  the  currents  are 
not  quantitatively  accurate.  However,  we  are  able  to  identify  those 

regions  of  the  magnetosphere  which  act  as  a  dynamo.  Moreover,  we  are  able 
to  at  least  partially  identify  the  driving  mechanisms. 

Method 

The  Birkeland  currents  were  derived  from  the  three-dimensional  MHD 
global  magnetosphere  simulation  model  [Brecht  et  al.,  1982].  The 

simulation  is  on  a  cartesian  mesh  with  a  resolution  of  two  earth  radii  near 
the  earth.  The  currents  are  computed  on  a  cubic  mesh  centered  on  the  earth 
with  the  faces  of  the  cube  located  at  about  five  earth  radii  geocentric 
altitude  as  shown  in  Figure  1.  There  are  36  mesh  points  on  each  face  of 
the  cube  so  that  76  mesh  points  map  into  each  of  the  north  and  the  south 
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polar  regions.  All  but  four  of  these  mesh  points  map  along  dipole  field 
lines  into  the  latitude  range  between  65  and  77  degrees. 

At  each  mesh  point  on  the  surface  of  the  cube  we  calculate  J  /Jb|  * 

(7  x  from  the  magnetic  field  intensity  B_  in  the  simulation 

results.  We  assume  that  within  the  cube,  the  flow  is  given  by  the 
electrostatic  approximation  with  infinite  field-aligned  conductivity.  Then 
the  quantity  J  / 1 B I  is  a  constant  along  flux  tubes  and  can  therefore  be 
easily  mapped  to  ionospheric  altitudes  along  the  assumed  dipole  field 
lines.  This  simplifying  assumption  is  probably  reasonable  during  most 
quasi  steady  flow  conditions.  The  magnitude  of  the  current  density  in  the 
ionosphere  is  then  obtained  by  multiplying  J(/|b|  by  the  ionospheric 
magnetic  field  intensity,  6  x  10“^  T.  Since  this  calculation  is  performed 
after  the  fact,  there  is  no  ionospheric  feedback  on  the  computed 
magnetosphere. 

For  the  results  shown  here  the  solar  wind  density  was  5  cm  J  and  the 
velocity  was  400  km  sec”^.  The  interplanetary  magnetic  field  was  5  nT 
southward.  The  magnetosphere  configuration  was  open  with  merging  between 
the  IMF  and  the  geomagnetic  field  occurring  on  the  bow  of  the 

magnetosphere.  Rapid  merging  was  not  occurring  in  the  magnetotail  since 
the  data  shown  was  selected  during  a  quiet  time  between  two  discrete 

reconnection  and  acceleration  events  in  the  tail. 

We  also  compute  an  example  of  the  plasma  convection  pattern  in  the 

polar  ionosphere.  The  electric  potential  is  derived  by  solution 

of  v»(I7i j>)  =  J(  ,  where  is  the  electric  potential  and  I  the  conductivity, 
was  chosen  to  be  a  uniform  scalar  5  mhos,  a  value  which  is  appropriate  to 
the  auroral  regions  but  too  large  in  most  of  the  polar  cap. 


Results 


Figure  2  shows  the  parallel  current  in  the  ionosphere  poleward  of  60 
degrees  latitude.  A  clear  Region  1  Birkeland  current  pattern  is 
apparent.  The  solid  contour  lines  show  the  upward  current  in  the 
afternoon-evening  sector  of  the  polar  region,  while  the  dashed  contours 
show  the  downward  current  on  the  morning  side.  No  clear  Region  2  current 
system  can  be  seen  primarily  owing  to  the  absence  of  data  points  below  65° 
latitude.  The  mid-day  currents  occur  primarily  at  latitudes  above  75°, 
while  at  dawn,  at  dusk,  and  at  night  the  currents  are  at  68°  to  72° 
latitude.  The  daytime  current?  are  concentrated  near  0900  and  1500  hours, 
the  evening  currents  between  1700  and  2300,  and  the  morning  currents 
between  0000  and  0700,  Current  density  magnitudes  range  generally  between 
0.1  and  0.2  x  10-^  amp  m-^. 

Figure  3  shows  the  polar  plasma  convection  consistent  with  the 
calculated  currents.  The  convection  potential  was  calculated  using  a 
uniform  ionospheric  conductivity  of  5  mho;  there  were  no  auroral 

conductivity  enhancements.  The  total  cross  polar  cap  potential  is  15  kV. 
The  convection  pattern  shows  the  typical  anti-solar  convection  poleward  of 
the  Region  1  current  system  and  sunward  flow  at  lower  latitudes. 

In  order  to  determine  the  source  regions  for  the  Region  1  current  in 
the  simulation  results,  we  have  traced  magnetic  field  lines  from  the  polar 
ionosphere  into  the  magnetosphere  and  solar  wind.  The  field  line  tracing 
results  show  a  clear  delineation  between  open  and  closed  geomagnetic 

field.  The  boundary  between  open  and  a  closed  field  lines  is  shown  in 

Figure  4  plotted  over  the  current  system  of  Figure  2. 

The  results  of  field  line  tracing  in  Figure  4  show  three  distinct 
source  regions  for  the  Region  1  Birkeland  currents.  The  dayside  currents, 


beCween  0800  and  1600  hours,  close  to  the  bow  magnetopause  and  to  the  polar 
magnetopause  on  open  polar  field  lines.  The  dawn  currents,  between  0400 
and  0700,  and  dusk  currents  between  1700  and  2000,  close  to  a  boundary 
layer  in  the  equatorial  plane  on  the  flanks  of  the  magnetosphere  and  are  on 
closed  field  lines.  The  nightside  currents,  between  2000  and  0400,  close 
along  the  boundary  between  the  plasma  shee't  and  the  tail  lobes  and  are  on 
open  field  lines.  The  characterization  of  an  individual  field  line  as  open 
or  closed  is  not  always  perfect  owing  to  limited  resolution  in  the 
simulation  and  to  numerical  error  in  the  field  line  tracing.  However, 
after  careful  study  of  a  number  of  different  cases  we  believe  that  the 
characterization  provided  above  is  accurate. 

Discussion  and  Conclusions 

In  this  discussion  we  would  like  to  address  a  number  of  concerns  and 
points  of  interest.  It  appears  that  the  Region  1  Birkeland  current  density 
is  about  one  order  of  magnitude  too  small  compared  to  observed  values 
[Iijima  and  Potemra,  1978),  whereas  the  integrated  current  densities  are 
about  a  factor  of  5  too  small.  The  differences  could  arise  from  a  number 
of  causes.  As  mentioned  previously,  the  numerical  resolution  (>  2  Rfi 
between  grid  points)  is  not  nearly  fine  enough  to  see  narrow  structures  in 
the  magnetosphere.  The  low  resolution  has  two  effects;  the  currents  tend 
to  be  spread  over  a  broader  area,  and  the  peak  magnitude  of  the  current  is 
reduced.  Unfortunately,  in  the  absence  of  higher  resolution  calculations 
it  is  difficult  to  assess  the  impact  of  resolution  on  these  results. 

Another  cause  of  the  reduced  Birkeland  current  magnitudes  is  related 


to  the  current  closure  mechanism  at  low  altitudes.  The  ionosphere  acts  as 
a  load  on  plasma  dynamics  in  the  magnetosphere.  As  the  ionospheric 


conductivity  is  increased  (decreased)  the  load  on  the  magnetosphere  is 
increased  (decreased)  and  the  driven  Birkeland  currents  increase  (decrease) 
in  magnitude.  The  magnetosphere-ionosphere  currents  are  regulated  by  a 
balance  between  the  strength  of  the  dynamo  and  the  size  of  the  load.  In 
the  simulation,  since  the  ionosphere  was  not  a  prescribed  boundary 
condition,  we  cannot  assess  the  size  of  the  load  and  therefore  cannot 

assess  the  Birkeland  current  magnitude  relative  to  observations.  In  the 
simulation  the  Birkeland  currents  close  at  low  altitude  via  a  loss  of 
resolution;  the  current  carrying  flux  tubes  cannot  be  distinguished  between 
widely  spread  grid  points.  Whether  this  closure  of  currents  in  some  sense 
resembles  the  ionosphere  conductivity  is  difficult  to  know.  However, 

current  closure  due  to  loss  of  numerical  resolutions  is  a  dissipative 
process.  As  such,  it  serves  as  a  load  on  the  dynamo  and  In  that  sense 
should  act  similarly  to  the  real  ionosphere. 

The  magnitude  of  the  convective  potential  across  the  polar  cap  of  15 
kV  is  also  smaller  than  the  average  observed  potential  of  about  40  kV  and 
should  probably  be  higher  for  a  5  gamma  southward  IMF  [Stern,  1979].  This 
difference  i  probably  due  to  both  the  smaller  values  of  the  sheet  current 
densities  and  the  simplified  model  for  the  ionospheric  conductivity  used  to 
generate  the  convection  pattern.  The  ionospheric  electric  fields  are  not 
fully  consistent  with  the  magnetospher ic  simulation  since  they  are 

calculated  after  the  fact  and  are  not  a  part  of  the  simulation. 

Nevertheless,  both  the  simulation  3irkeland  currents  and  convection 
potentials  bear  a  strong  resemblance  to  those  actually  observed. 

We  now  turn  our  attention  to  the  current  source  regions.  The  boundary 
between  open  and  closed  geomagnetic  field  lines  clearly  divides  the  Region 
l  currents  into  three  distinct  systems.  At  dawn  and  at  dusk  the  currents 


flow  on  closed  field  lines,  while  during  the  day  and  at  night  the  currents 
flow  either  along  the  boundary  or  on  open  field  lines.  These  systems  will 
now  be  discussed  individually  in  greater  detail. 

The  geomagnetic  field  lines  which  carry  the  dawn  and  dusk  currents  are 
highly  distorted  in  the  anti-solar  direction.  Field  lines  poleward  of  the 
current  sheet  and  penetrating  the  ionosphere  at  0700  and  1700  hours  cross 
the  equatorial  plane  about  2  earth  radii  in  the  anti-solar  direction  while 
those  near  0400  and  2000  hours  extend  more  than  40  earth  radii  tailward. 
Field  lines  equatorward  of  the  sheet  are  substantially  less  distorted.  All 
these  current  carrying  field  lines  cross  the  equatorial  plane  in  a  boundary 
layer  near  the  magnetopause.  The  reader  might  expect  that  a  low  latitude 
boundary  layer  current  would  extend  further  toward  mid-day  than  0700  and 
1400  hours.  In  reality  they  probably  do  and  the  termination  in  the 
simulation  results  is  associated  with  numerical  resolution  problems  along 
the  dayside  magnetopause. 

The  current  systems  on  closed  field  lines  are  driven  by  a  momentum 
coupling  process  from  the  magnetosheath,  across  the  magnetopause,  and  into 
a  low  latitude  magnetosphere  boundary  layer.  In  the  simulation  code  this 
coupling  is  due  to  both  numerical  viscosity  and  an  artificial  viscosity  in 
the  momentum  equation  [Brecht  et  al.,  1982].  At  the  present,  however,  we 
are  unable  to  make  a  quantitative  comparison  between  the  simulation 
viscosity  and  an  actual  physical  viscosity  in  the  magnetosphere.  The 
effective  viscosity  in  the  simulation  is  likely  to  be  substantially  larger 
than  the  viscosity  which  actually  pertains  to  the  low  latitude  boundary 
layer.  We  intend  to  address  this  question  in  a  more  quantitative  manner  in 
the  future. 
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The  dayside  and  nightside  currents  have  a  different  driving 
mechanism.  It  is  tempting  to  associate  the  currents  on  open  field  lines  to 
a  process  which  is  driven  by  magnetic  merging.  This  appears  not  to  be  the 
case.  During  the  simulations  we  have  observed  sporadic  rapid  merging  and 
acceleration  of  plasma  in  the  magnetotail.  The  merging  and  acceleration 
processes  occur  with  virtually  no  disturbance  of  the  Region  1  current 
systems.  We  therefore  conclude  that  the  current  systems  presented  here  are 
not  associated  with  magnetic  merging  processes. 

The  dayside  Birkeland  currents  on  open  field  lines  appear  to  cross  the 
nose  magnetopause  and  extract  momentum  from  the  anti-solar  magnetosheath 
flow.  By  loading  the  currents  in  the  magnetosheath  which  reaccelerate 
plasmas  away  from  the  bow  stagnation  point  the  Birkeland  currents  and  the 
ionosphere  draw  energy  directly  from  the  solar  wind.  These  dayside  open 
field  line  currents  appear  to  occur  in  the  same  region  as  the  cusp  currents 
of  Iijima  and  Potemra  [1976b]  and  D'Angelo  [1980],  Unfortunately,  the 
simulation  currents  have  opposite  polarity  to  the  cusp  systems  which  were 
observed.  We  are  unable  to  explain  this  discrepancy. 

The  simulation  nightside  currents  occur  in  the  same  region  as  those 
observed  in  the  magnetosphere  [Aubrey  et  al.,  1972;  Fairfield,  1973; 
Sugiura,  1975;  Frank  et  al.,  1981;  Kelley  et  al.,  1981].  The  mechanism 
responsible  for  driving  these  nightside  currents  has  not  been  identified  at 
this  time.  The  analysis  of  the  source  region  is  very  difficult  owing  to 
the  weak  fields,  low  plasma  densities,  apparently  unorganized  flow,  and  low 
numerical  resolution  in  this  region  of  the  simulation.  It  is  possible  that 
the  mechanism  of  Rostoker  and  Bostrora  [1976]  may  be  responsible  for  these 
currents  but  at  present  we  cannot  be  sure.  Study  is  continuing  to  identify 
the  source  of  these  currents. 


In  conclusion,  we  would  like  Co  summarize  Che  resulcs.  Numerical  MHD 
siraulacions  of  Che  magnecosphere  produce  3irkeland  currenc  syscems  which 
are  in  reasonable  agreemenc  wich  observaCions .  Here  we  have  presenced 
resulcs  for  Che  Region  1  currencs.  These  currenCs  are  divided  by  Che  open- 
closed  field  line  boundary  inco  Chree  separace  sycems;  a  low  lacicude 
boundary  layer  syscem  along  Che  flanks  of  Che  magnecosphere,  a  dayside  open 
field  line  syscem  which  has  ics  source  in  Che  magneCosheaCh ,  and  a 
nighcside  open  field  line  syscem  on  Che  boundary  beCween  Che  Call  lobes  and 
che  plasma  sheec.  We  have  also  discussed  possible  source  mechanisms 
responsible  for  Chese  syscems. 

The  resulcs  presenced  here  are  preliminary.  Improvemencs  are  being 
made  in  Che  simulacion  codes  with  regard  to  resolution  and  in  order  to 
simulate  a  realistic  ionosphere  boundary  condition.  A  factor  of  two 
resolutions  has  been  obtained  by  refinement  of  the  numerical  mesh,  and  an 
additional  factor  of  two  has  been  achieved  by  upgrading  the  computational 
algorithms.  The  new  ionospheric  boundary  conditions  allow  ionospheric 
fields  and  flows  to  feed  back  on  the  magnetospheric  dynamics.  We  expect 
that  the  continued  work  will  lead  to  improved  results  in  the  near  future. 
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CUBIC  MESH  FOR  Jn 


A  diagram  of  the  cubic  mesh  used  to  calculate  parallel  currents 
from  the  global  simulation  results.  The  earth  is  located  at 
the  center  and  each  face  of  the  cube  is  about  5  earth  radii 
from  the  earth.  The  currents  are  mapped  from  the  mesh  points 
on  the  cube  along  dipolar  field  lines  to  the  polar  ionosphere. 


Solid/dashed  contours  are  positive/ 
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P.O.  BOX  969 
LIVERMORE,  CA  94550 

0 1 C  Y  ATTN  DOC  CON  FOR  8.  MURPHEY 

01 C  Y  ATTN  DOC  CON  FOR  T.  COOK 
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OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

01CY  ATTN  DOC  CON  DR.  YO  SONG 


OTHER  GOVERNMENT 


BOEING  COMPANY,  THE 
P.  0.  8  0  X  3707 
SEATTLE,  WA  98124 

01 CY  ATTN  G.  KEISTER 
0 1 C  Y  ATTN  D.  MURRAY 
0 1  C  Y  ATTN  G.  HALL 
0 1 C  Y  ATTN  J.  KENNEY 


INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  S  INFO 
ADMIN 

BOULDER,  CO  80303 

0 1 C  Y  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01 CY  ATTN  U.  UTLAUT 
0 1  C  Y  ATTN  D.  CR0M3IE 
0 1 C Y  ATTN  L.  BERRY 

NATIONAL  OCEANIC  S  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
0 1 C  Y  ATTN  R.  GRUBB 
01  CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01 CY  ATTN  I.  GARFUNKEL 
0 1 C  Y  ATTN  T.  SALMI 
01 CY  ATTN  V.  JOSEPHSON 
01 CY  ATTN  S.  80WER 
01 CY  ATTN  D.  OLSEN 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 

555  TECHNOLOGY  SQUARE 

CAMBRIDGE,  MA  02139 
01CY  ATTN  D.B.  COX 
01 CY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 

LINTHICUM  ROAD 

CLARKSBURG,  MD  20734 
01  CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

ITHACA,  NY  14850 

01  CY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
01  CY  ATTN  H.  LOGSTON 
01 CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wi  l s h i re  B  L  vd . 

Santa  Monica,  Calif  90401 
01 CY  ATTN  C.B.  GABBARD 
01 CY  ATTN  R.  LELEVIER 


ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

0 1 C  Y  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1 901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01 CY  ATTN  L.  SLOAN 
01 CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 

01 CY  ATTN  J.  WORKMAN 
0 1 C  Y  ATTN  C.  PRETTIE 
01 C  Y  ATTN  S.  BRECHT 


ESL ,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
01 CY  ATTN  J.  ROBERTS 
01 CY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
01 CY  ATTN  M.H.  BORTNER 
SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 

01 CY  ATTN  F.  RE  I  BE  RT 


y.O'.V.V/ 


GENERAL  ELECTRIC  TECH  SERVICES 
CO . ,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

0 1  C  Y  ATTN  G.  HILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIR3ANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01 CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01 C  Y  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTSRN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01 C  Y  ATTN  DICK  STEINHOF 

HSS ,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01 CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01 CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  V A  22311 
0 1 C  Y  ATTN  J.M.  A  E I N 
OIC*  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  UOLFARD 
01CY  ATTN  JOEL  3ENGST0M 

INTL  TEL  8  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

0 1 C  Y  ATTN  TECHNICAL  LIBRARY 

J  AY  COR 

11011  TORREYANA  ROAD 

P.O.  BOX  35154 

SAN  DIEGO,  CA  92138 

01CY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20310 

0 1  C  Y  ATTN  DOCUMENT  LIBRARIAN 
0 1  C  Y  ATTN  THOMAS  POTEMRA 
0 1 C  Y  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01 C  Y  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
0 1 C  Y  ATTN  DAS  I  AC 
01 CY  ATTN  WARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
0 1 C  Y  ATTN  8.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
0 1  C  Y  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01  CY  ATTN  DEPT  60-12 
01  CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-12 
01  CY  ATTN  W.L.  IMHOF  DEPT  52-12 
0 1 C  Y  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01 CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32S05 

01CY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 C  Y  ATTN  DAVID  M.  TOWLE 
01 CY  ATTN  L.  LOUGHLIN 
0 1 C  Y  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  50LSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 


01  CY 

ATTN 

N .  HARRIS 

01  CY 

ATTN 

J.  MOULE 

01  CY 

ATTN 

GEORGE  MROZ 

01  CY 

ATTN 

W.  OLSON 

01  CY 

ATTN 

R.W.  HALPRIN 

01  CY 

ATTN 

TECHNICAL 

LIBRARY  SERVICES 

MISSION 

RESS 

ARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 

01  CY 

ATTN 

P.  FISCHER 

01  CY 

ATTN 

W.F.  CREVIER 

01  CY 

ATTN 

STEVEN  L.  GUTSCHE 

01  CY 

ATTN 

R.  BOGUSCH 

01  CY 

ATTN 

R.  HENDRICK 

01  CY 

ATTN 

RAl-PH  kilb 

01  CY 

ATTN 

DAVE  SOWLE 

01  CY 

ATTN 

F.  FAJEN 

01  CY 

ATTN 

M.  SCHEIBE 

01  CY 

ATTN 

CONRAD  L.  L0N3MIRE 

01  CY 

ATTN 

B.  WHITE 

01  CY 

ATTN 

R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NEW  MEXICO  37106 
01  CY  R.  STELLINGWEP.F 
0 1 C  Y  M.  ALME 
01 C  Y  L.  WRIGHT 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  L A3 
31 S  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMSTRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

0 1  C  Y  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  93009 

01CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  S  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


01  CY 

ATTN 

FORREST  GILMORE 

01  CY 

ATTN 

WILLIAM  B.  WRIGHT 

01  CY 

ATTN 

WILLIAM  J.  KARZAS 

01  CY 

ATTN 

H.  ORY 

01  CY 

ATTN 

C.  MACDONALD 

01  CY 

ATTN 

R.  TURCO 

01  CY 

ATTN 

L.  DeRAND 

0 1  C  Y 

ATTN 

W.  TSAI 

MITRE  CORPORATION,  THE 
P.O.  BOX  203 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  JOHN  MORGANSTSRN 
01 CY  ATTN  G.  HARDING 
01 CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  V A  22101 
01 C  Y  ATTN  W.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 C  Y  ATTN  E.C.  FIELD,  JR. 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
01 CY  ATTN  ED  BED  ROZ IAN 

RAYTHEON  CO. 

528  80ST0N  POST  ROAD 
SUDBURY,  MA  01776 

01 CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01 C  Y  ATTN  VINCE  TRAPANI 


25 


<r^  •.  «.  ««  _  •  .  «  »  «  ■■  •  k  «•  -  .  *»  *■  ■•  *  •  *  *  *  *  * 


VvT 


,  i ». 


% , 

? 

!%. 

$ 

* 

i  • 

S' 

I 

<r 


SCIENCc  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

0 1 C  Y  ATTN  LEWIS  M.  LINSON 
0 1 C  Y  ATTN  DANIEL  A.  HAMLIN 
0 1 C  Y  ATTN  E.  FRIEMAN 
0 1 C Y  ATTN  E . A  .  STRAKER 
01 CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  G00DRI0GE  DR. 

MCLEAN,  VA  22102 
01 CY  J.  COCKAYNE 
0 1 C  Y  E.  HYMAN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 


MENLO 

PARK, 

CA  94025 

01  CY 

ATTN 

J.  CASPER 

01  CY 

ATTN 

DONALD  NEILSON 

01  CY 

ATTN 

ALAN  BURNS 

01  CY 

ATTN 

G.  SMITH 

01  CY 

ATTN 

R.  TSUNODA 

01  CY 

ATTN 

DAVIi'  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHE3NUT 

Cl  CY 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

01  CY 

ATTN 

J.  VICKREY 

01  CY 

ATTN 

RAY  L.  LEADA3RAND 

01  CY 

ATTN 

G.  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

01  CY 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V.  GONZALES 

01  CY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01 CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01 CY  ATTN  JOHN  ISE,  JR. 

01 CY  ATTN  JOEL  GAR9ARIN0 

TRW  DEFENSE  8  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90273 
01 CY  ATTN  R.  K.  PLE2UCH 
0 1 C  Y  ATTN  S.  ALTSCHULER 
01 CY  ATTN  D.  DEE 
01 CY  ATTN  0/  STOCKWELL 
SNTF/1575 


V  I S  I  D  Y  N  E 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01 CY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
01 CY  ATTN  C .  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01 C  Y  ATTN:  N.  ZABUSKY 
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IONOSPHERIC  MODELING  DISTRI3UTI0N  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
Code  4100 

Dr.  H.  GURSKY  -  CODE  4100 
Dr.  P.  GOODMAN  -  CODE  4180 

A. F .  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01731 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  W.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERS 
Dr.  F.J.  RICH 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  14350 
DR.  W.E.  SWARTZ 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 

PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DUB IN 

DR.  N.  MAYNARD  -  CODE  696 
COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 

ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L. A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D. A.  HAMLIN 
DR.  E.  FRIEMAN 
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STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
OR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
DR.  L.E.  LEE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
DR.  M.  PONGRATZ 
DR.  D.  SIMONS 
DR.  G.  8ARASCH 
DR.  L.  DUNCAN 
DR.  P.  BERNHARDT 
DR.  S.P.  GARY 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

J-OHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20310 
DR.  R.  GREENWALD 
DR.  C.  MENG 


UNIVERSITY  OF  TEXAS 
AT  DALLAS 

CENTER  FOR  RESEARCH  SCIENCES 
P.O.  BOX  688 
RICHARDSON,  TX  75080 
DR.  R.  HEELIS 
DR.  W.  HANSON 
DR.  J.P.  McCLURE 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 
DR.  J.  ST. -MAURICE 

PHYSICAL  RESEARCH  LABORATORY 
PLASMA  PHYSICS  PROGRAMME 
AHMEDABAD  380  009 
INDIA 

P.J.  PATHAK,  LIBRARIAN 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERGY  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
JHAN  VARYAN  HELLMAN, 
REFERENCE  LIBRARIAN 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 
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